Formation time scaling and hadronization in cold nuclear matter 
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I propose a scaling analysis of the hadron multiplicity ratio measured in Deep Inelastic Scattering 
on nuclear targets as a tool to distinguish energy loss and nuclear absorption effects on hadron sup- 
pression in cold nuclear matter. The proposed scaling variable is a function of the hadron fractional 
energy and of the virtual photon energy. Its functional form, which depends on a parameter A, 
can be fixed by general theoretical considerations and encompasses both energy loss and absorption 
models. The parameter A is fitted to HERMES experimental data and shown to favor prehadron 
nuclear absorption as leading mechanism for hadron suppression as opposed to quark energy loss. 
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I. INTRODUCTION 

In Deep Inelastic Scattering on nuclear targets (nDIS) 
one observes a suppression of hadron production [1-7] 
analogous to hadron quenching in heavy-ion collision at 
the Relativistic Heavy-Ion Collider (RHIC) [8]. 

The cleanest environment to address nuclear modifi- 
cations of hadron production is nuclear DIS: it allows 
to experimentally control many kinematic variables; the 
nuclear medium, i.e., the nucleus itself, is well known; 
the multiplicity in the final state is low. Moreover, the 
nucleons act as femtometer-scale detectors allowing to 
experimentally study the propagation of a parton in this 
"cold nuclear matter" , and its space-time evolution into 
the observed hadron. In the case of heavy ion collisions, 
one wants to use hadron suppression as a tool to extract 
the properties of the hot and dense system created in 
the collision, also called "hot nuclear matter" . If, for 
example, the parton's color were neutralized on much 
larger scales than the nuclear radius, hadron suppression 
would be attributed to parton energy loss [9]. Analysis 
of midrapidity hadron production at RHIC in the energy 
loss framework leads to a medium temperature T w 400 
MeV, well in excess of the critical temperature Tc ~ 170 
MeV for the transition into a deconfined Quark-Gluon 
Plasma [10, 11]. If, on the contrary, color neutralization 
started on the nuclear radius scale or before, one should 
also account for the interactions of the medium with the 
prehadron, the color neutral precursor of the hadron [12]. 
This would lead to a different, presumably lower, value 
of the medium temperature. Knowing precisely how the 
struck quark propagates in cold nuclear matter - most 
importantly, whether it starts hadronizing inside or out- 
side the nuclear medium - is essential for correctly using 
hadron quenching as a signature of the production of a 
Quark-Gluon Plasma at RHIC. 

Experimental data on hadron production in nDIS are 
usually presented in terms of the multiplicity ratio [1-6] 



i.e., the single hadron multiplicity per DIS event on a 
target of mass number A normalized to the multiplicity 
on a deuterium target, as a function of the virtual pho- 
ton energy ly and of Zh = p ■ Ph/p • <?, with p the target 
4- momentum divided by A, ph the hadron 4-momentum 
and q the virtual photon 4-momentum. In the target rest 
frame Zh = Eh/v is the hadron fractional energy with re- 
spect to the virtual photon energy. The double ratios 
in (1) cancel to a large extent initial state effects like 
the modifications of parton distribution functions due to 
shadowing and EMC effects, exposing the nuclear mod- 
ifications of the fragmentation process. If no nuclear ef- 
fects modified the fragmentation process, we would ex- 
pect Rm ~ 1. In fact, what is experimentally observed 
[1-4] is a suppression of pions, kaons and antiprotons in 
the Zh = 0.1 — 1 and v = 1 — 100 GeV range. Protons are 
enhanced at z/i < 0.4 ("proton anomaly") and suppressed 
above. Both quenching and enhancement increase with 
A. 

Despite a lot of experimental and theoretical efforts, 
the leading physical mechanism for hadron quenching in 
nDIS has not yet been unambiguously established. In 
particular, as shown in [13, 14], the observed approximate 
A^/^ scaling of the experimental data cannot distinguish 
models based on nuclear absorption [13-18] from models 
based on parton energy loss [11, 19], as is often assumed. 
Indeed, single hadron suppression in nDIS obeys a A^^^ 
law (broken at A > 80) in both energy loss and absorp- 
tion models [13]. Even the more refined analysis in terms 
of Rm — cA" fits proposed in [14] cannot clearly distin- 
guish the 2 classes of models. 

In this paper, I propose a scaling analysis of Rm as 
a tool to disentangle parton energy loss and nuclear ab- 
sorption effects on hadron production in nDIS. More in 
detail, I conjecture that Rm should not depend on z^ 
and v separately but should depend on a combination of 
them: 
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Rm = RM[T{zh,iy)] , 
where the scaling variable r is defined as 

T = C Z^{1 - Zh)iy . 



(2) 



(3) 



The scaling exponent A is introduced as a way of approx- 
imating and summarizing the scahng behavior of experi- 
mental data and theoretical models. It will be separately 
obtained by a best fit analysis of data and theoretical 
computations, see Section III. The proportionality con- 
stant C cannot be determined by the fit. A possible 
scaling of Rm with Q^ is not considered in this analysis 
because of its model dependence; moreover, in the HER- 
MES data considered in this paper, the dependence of 
the average (Q^) on Zh and v is very mild, implying very 
small effects on the scaling of Rm- 

As discussed in Section II, the proposed functional 
form of T is flexible enough to encompass both absorption 
models and energy loss models. The 2 classes of models 
are distinguished by the value of the scaling exponent: 
a positive A ^ is characteristic of absorption models, 
while a negative A < is characteristic of energy loss 
models. Thus, the exponent A obtained in the proposed 
model-independent scaling analysis of experimental data 
can identify the leading mechanism for hadron suppres- 
sion in nDIS. 



II. SCALING OF Rm 

The idea that the hadron multiplicity ratio Rm should 
scale with the variable r introduced in Eq. (3) is quite 
natural in the context of hadron absorption models [14- 
18]. In these models the struck quark neutralizes its color 
on a relative short time scale. The ensuing color neu- 
tral state, called a prehadron, later on collapses on the 
wave function of the observed hadron. Hadron suppres- 
sion is then mainly attributed to prehadron-nuclcons in- 
teractions, whose magnitude depends on the in-medium 
prehadron path length, which depends solely on the pre- 
hadron formation time i*. 

Estimates of the prehadron formation time can be 
obtained in the framework of the Lund string model 
[14, 15, 18], where the prehadrons are identified with 
each of the fragments of the color string. Alternatively, in 
the pQCD inspired color dipole model for leading hadron 
suppression of Ref. [16] , the hardest gluon radiated off the 
struck quark splits into a quark-antiquark pair; the an- 
tiquark then recombines with the struck quark into the 
leading prehadron. In both cases the prehadron forma- 
tion time has a simple general form: 



t* = g{zh){l - Zh)- 

K 



(4) 



where ^(z^)— »0 as z/j— >0, and k is a constant that sets 
the time scale of hadronization. In the Lund model k « 1 
GeV/fm is given by the string tension; in the dipole 
model K = Q2. At HERMES, in both models, the pre- 
hadron formation time is <* < 5 fm, which is smaller 
than the nuclear radius. On the contrary, hadrons are 
typically produced at the periphery or outside the target 
nucleus and their absorption does not contribute much 
to Rm- The physical origin of t, is transparent. The 



factor v can be understood as a Lorentz boost factor. 
At large Zh the hadron carries away most of the struck 
quark energy. The color string remainder has only an 
energy (1 — Zh)v left, so that it cannot stretch farther 
off (in pQCD terms, the colored struck quark has a little 
energy to radiate into gluons, hence it must neutralize its 
color in a short time). At small z^^O the prehadron for- 
mation time should go to 0, as well, as explicitly shown 
in Lund model computations [14, 15, 18]. This follows 
from the fact that we are discussing semi-inclusive hadron 
measurements. At small Zh the observed hadron carries 
away a small fraction of the struck quark energy. The 
rest of the energy will most probably be used for the cre- 
ation of other low energy prehadrons, because the string 
fragmentation function is steeply falling with Zh- On av- 
erage, the observed prehadron will be produced close to 
the interaction point. 

Summarizing the above discussion, Rm in absorption 
models should depend only on t* = t*(z/i, i^) and not on 
Zfi and ly separately. A good approximation to t* is the 
scaling variable r of Eq. (3), where the scaling exponent 
A depends on the chosen absorption model. A rough 
estimate of the scaling exponent gives A « 1. A more 
precise value can be obtained by fitting Eqs. (2)-(3) to 
the theory model results [13, 14, 16, 17] for Rm- The fit 
procedure, explained in detail in the next section, results 
in 0.5 < A < 1.2 for absorption models. 

In energy loss models [11, 19] the hadron formation 
time is assumed to be much larger than the nuclear ra- 
dius, and the hadronization process is assumed to happen 
entirely outside the target nucleus [25] . The quark travels 
through the nucleus and experiences multiple scatterings 
and induced gluon bremsstrahlung. Hence, it starts the 
hadronization process with a reduced energy v — e where 
e is the energy of the radiated gluons 

In Ref. [19], extended in [13] to include finite medium 
size corrections, the reduced quark energy at the time of 
hadronization is translated into a shift of Zh in the vac- 
uum fragmentation function D [20] . The medium modi- 
fied FF is then computed as 
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where the dependence of the vacuum FF on the hard scale 
scale Q^ of the process is understood, and the quenching 
weight 'P(e) is the probability distribution of an energy 
loss e computed in the Baier-Dokshitzer-Muellcr-Schiff 
formalism [21]. Note the upper limit of integration in 
Eq. (5) imposed by energy conservation. For the pur- 
pose of discussing the scaling properties of Rm, we can 
work in the soft gluon approximation, and neglect finite 
quark energy corrections, which would introduce an ad- 
ditional V dependence in the quenching weight [22]. If we 
further neglect energy loss fluctuations, we can approxi- 
mate Rm ~ DA{zh)/D(zh) and obtain 
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where the average energy loss (e) = 

J^'-^'^>deene)/J^'-'^^''dene) = /[(I - z,)H is 
a function of the energy (1 — z^)!^ not carried away by 
the observed hadron. Next, we can approximate the FF 
using the parametrization of Ref. [23] at Q^ = 2 GeV^: 
D{zh) — Cz'jl{\ — ZhY , where for pions a ~ —1, (3 « 1.5 
and the constant C will cancel in the multiplicity ratio. 
Finally, 
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which shows an approximate scaling with (1 — Zh)^. 

In Ref. [11] the medium modifications of the fragmen- 
tation functions are computed from twist-4 contributions 
to the leading order cross-section, including diagrams 
with one elastic quark-nucleus scattering and one radi- 
ated gluon. Both the struck quark and the radiated gluon 
are allowed to fragment according to vacuum FF. The ob- 
tained modified FF, D, can be well approximated by the 
numerator in Eq. (6) with t/v ~ 0.6(zg}, where {zg) is the 
average fractional energy of the radiated gluon [11, 24]: 
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Here, Jq^gg is the quark-gluon splitting function, C{Q^) 
is the strength of parton-parton correlations in the nu- 
cleus, rriN the nucleon mass, Ra the nuclear radius, and 
fc is a shorthand for the quantities independent of Zh and 
p. fg is s function of 1 — z^ because of the upper limit 
of integration on Zg imposed by energy conservation. In 
the HERMES regime, /g(l - Zh) oc (1 - Zh)°-'^. Approxi- 
mating Rm and the modified FF as before we have: 
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From Eqs. (7) and (9) a scaling of Ri\i with (1 — Zh) 
is evident, which implies A = in Eq. (3). However, 
it is not immediate to see the role played by v. To 
establish it, let's introduce an effective scaling variable 
r' = Cz^{l — Zh)v^ , with /i an effective parameter de- 
scribing the scaling of Rm with respect to v in energy 
loss models. The value of /i can be determined by a fit 
of the full computations in Refs. [11, 13] as follows. For 
any given /i, we fit the theoretical Rm = Rm{t') and 
determine A — X{p) by x^ minimization as described in 
Section III. A scaling of Rm with r' (i.e., x^/d.o.f. < 1) 
is found for 0.2 < /x < 1.8, with the best-fit A decreasing 
as fi increases. As also expected on theoretical grounds. 
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FIG. 1: An example of the fit procedure described in Section 
3, applied to tt"*" production on a Kr target at HERMES [3]. 
Upper panel: x^ ^^ ^ function of A. Lower panel: Rm{t) with 
r computed at Abcst = 0.34. Experimental statistical errors 
are of the same size as the plotted points. 



in this range of /i values one finds Abost ^ 0, which dis- 
tinguishes it from the positive A expected in absorption 
models. 

In conclusion, for the sake of comparing energy loss 
models with absorption models, where /i = 1 is theoreti- 
cally justified, I will fix /i = 1 and analyze experimental 
data and theory models in terms of the scaling variable 
T proposed in Eq. (3). 

Finally, a note on the limitations of the proposed scal- 
ing. In absorption models, the proposed scaling might 
be broken by the dependence of the prehadron cross sec- 
tion cr* on the photon virtuality Q^ and the prehadron 
energy _E* w z^iy, and by a possible dependence of i* 
on Q^ [16]. At present, all these effects cannot be cal- 
culated from first principles, and are to a good extent 
model dependent. The Q^ dependence of a^: and t.^, is 
not a concern for the analysis of HERMES data because 
the range of (Q^) in their Zh- and I'-bins is rather small, 
implying a small scale breaking. In future experimental 
analyses it will be important either to measure Rm at 
fixed Q^ or to ensure that (Q^) stays approximately con- 
stant in all bins. Since the prehadron must evolve into 
the observed hadron on a relatively short timescale, one 
may expect that cr* ex ah{Eh) when averaging ct* along 
the prehadron path [14, 15]. Since au has a mild depen- 
dence on the hadron energy Eh in the HERMES kine- 
matic regime, only minor deviations from scaling are ex- 
pected. In Ref. [16], ct* is computed in the pQCD dipole 
model and is explicitly Eh and Q^ dependent. The good 
fit of A obtained in Section IV for this model is an a pos- 
teriori indication that scale breaking effects are small. 
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FIG. 2: The scaling exponent Abost extracted from HERMES 
data on charged and identified hadrons at Eiab = 27 GeV 
[2-4] (only statistical errors included in the fit). Error bars 
correspond to 1 standard deviation. Bottom panel; y^ per 
degree of freedom. 



In energy loss models, scale breaking may arise due to 
fluctuations in energy loss, especially near the kinematic 
limit. The relatively large error bars on A found in the 
fits of Section IV show that their effect is not fully negli- 
gible [26]. However, this docs not spoil the discriminative 
power of A, which yields A < for energy loss models but 
A ^ for absorption models. 



III. FIT PROCEDURE 

The HERMES experiment measures Rm binned in Zh 
and integrated over v and Q^ {"zh distributions") or 
binned in u and integrated over Zh and Q^ ("i^ distri- 
butions" ) . The scaling of experimental data with respect 
to the variable r defined in Eq. (3) and the scaling expo- 
nent A can be determined by a fit to the data as follows. 

(1) Fix A. 

(2) For each z^ bin in z/j-distributions compute r = 
T{zh, {y{zh))) andRMir) = Rpiizh), where {v{zh)) 
is the average measured v in the considered Zh-hin. 
Likewise for each v bin in (^-distributions compute 
T = T{{zh{iy)), v) and Rm{t) = Rm{i^)- 

(3) Fit a function (/)(t) to the pairs {(t, Rm)} obtained 
at step 2, and compute x^ = X^{^)- The choice of 
(j) is discussed below. 

(4) Determine the best-fit exponent Abost by minimiza- 
tion of x^(A). 

(5) If x^(Abost) ^ 1 per degree of freedom, we say that 
the analyzed data set scales with respect to t and 
is characterized by a scaling exponent Abcst- 
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FIG. 3: The scaling exponent Abcst for HERMES data at iSiab 
=12 GeV on a Kr target [5]. The N target results are not 
shown because of their very large error bars. Bottom panel: 
X^ per degree of freedom. 



An example of this procedure and the corresponding 
Rm{t) computed at A = Abcst is illustrated in Fig. 1. The 
fit to theoretical computations is done in the same way 
as the fit to HERMES data, by considering the computed 
Rm at the central value of each Zh, and v experimental 
bin. Theoretical errors are estimated as 6% of 1 — Rm 
for the models of Refs. [11, 13, 14, 17], which need to fit 
1 parameter to Rm data [14], and 10% for the model of 
Ref. [16]. 

The fit results discussed below have been obtained us- 
ing as fit function (/i(t) a polynomial of 4**^ degree in r. 
The results of the fit have been cross-checked by using a 
rational function of second order constrained to tend to 
1 as T-^oo, which has 5 free parameters as the default 
polynomial. A second cross-check was obtained by ad- 
ditionally constraining the rational function to have null 
derivative at r = in order to avoid singularities. Fi- 
nally, unless otherwise explicitly stated, for this scaling 
analysis I considered only data points satisfying the fol- 
lowing cuts, (i) Zh > 0.2 {{zfi) > 0.2 for t^-distributions), 
to avoid the target fragmentation region and feed-down 
of hadrons from higher- z;i, for which the conjectured scal- 
ing is not valid, and to avoid large corrections due to the 
detector geometric acceptance [17]. A cut at z^ > 0.3 
or 0.4 might be preferable from this point of view, but 
excessively reduces the available number of data points. 
(ii) Zh < 0.9, to avoid diffractive hadron production and 
quasi-elastic lepton-nucleus scatterings, (iii) v > 7 GeV, 
for consistency between the analysis of the N, Kr and 
Ne target data sets. For each target and hadron flavor 
15 data points survive these cuts. I explicitly checked 
the stability of A against small variations of the upper Zh 
cut. The stability of A against variations of the lower Zh 
cut is difficult to establish because of a rapidly shrinking 
number of data points with increasing {zh)min- 



IV. RESULTS 

The scaling exponents Abcst extracted from HERMES 
data at Eiab ^ 27 GeV [2-4] and 12 GeV [5] for differ- 
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FIG. 4: Comparison of the scaling exponent for vr and h 
from HERMES data at Eub = 27 GeV [2-4] and from theory 
models. Error bars correspond to 1 standard deviation. En- 
ergy loss models (blue points on-line): AA [13], WW [11]. Ab- 
sorption models (red points): AGMP (pure absorption with- 
out Q'^-rescaling) [13, 14], Col. Dip. [16]. The Giessen model 
[17] embeds nuclear absorption in a full Monte Carlo simu- 



lation of the nDIS event, 
freedom. 



Bottom panels: x P^r degree of 



ent hadron flavors produced on N, Ne and Kr targets are 
shown in Fig. 2. In all cases x^/d.o.f. < 1.6, which proves 
that Rm scales with r. The central result of this paper is 
that pion data exhibit a clear Abost ^ 0. As discussed in 
Section II, this result contradicts the assumption used in 
energy loss models that the quark is long-lived. In other 
words, it indicates the dominance of the prehadron ab- 
sorption mechanism as opposed to the energy loss mech- 
anism. This conclusion is confirmed by the comparison 
to theory models shown in Fig. 4. Therefore, the scaling 
variable r can be interpreted as a measure of the pre- 
hadron formation time. 

Unidentified charged hadrons {h ) have a positive 
Abcst on N target, but Abcst ~ on Kr target. This ap- 
parently contradictory result can be explained in terms 
of the proton contribution to the h sample. Proton pro- 
duction shows an anomalous enhancement of Rm above 
1 when Zh ^ 0.4, which cannot be explained in terms 
of either parton energy loss or prehadron absorption, 
and is not yet fully understood theoretically. The pro- 
ton anomaly explains the negative value of its scaling 
exponent Abcst, which in turn drives the h value of 
Abcst towards for heavy targets. Indeed, by cutting 
the h^ sample on Kr target at z^ > 0.5 one obtains a 
reduced x^/d.o.f ~ 1.18 and Abcst — 0.34 ± 0.13, com- 
patible with the TT exponents and the nuclear absorp- 
tion mechanism. A further confirmation of the role of 
proton anomaly in reducing the scaling exponent comes 
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FIG. 5: Scaling exponent and x^ per degree of freedom for 
"1ow-t" and "high-r" data sets defined in Section IV, com- 
pared to the full data set for pion production on Ne targets, 
corresponding to "medium-r". Black disks: -k^ from prelim- 
inary HERMES data [4]. Red circles: 7r+ from the AGMP 
absorption model [13, 14]. Blue triangles: tt from the en- 
ergy loss model of Ref. [13, 19]. (Color on line). 



from preliminary data on Kr at Eiab = 12 GeV [5], which 
yield Abcst = 0.02 ± 0.09 for h+ but Abcst = 0.34 ±0.11 
for h^ , see Fig. 3. From kaons and antiprotons data it 
is difficult to draw any conclusion because of the large 
error bars. 

An interesting cross-check of the interpretation of t 
as the prehadron formation time can be obtained by di- 
viding the full data set in two subsets with low- and 
high-r. In the high-r data set the prehadron has a 
shorter in-medium path length, and the quark has a 
longer in-medium quark path length, than in the low- 
r data set. Then, one expects a smaller contribution 
of prehadron absorption and a larger contribution from 
partonic interactions and energy loss, hence a smaller A 
[27]. Since r monotonically decreases with Zh in HER- 
MES z^-distributions, and monotonically increases with 
1/ in z^-distributions, we can define the 2 subsets by the 
following cuts on Zh and v. 

• low-r: 0.5 < z < 0.9 and 7 GeV < t/ < 13 GeV 

• high-r: 0.2 < z < 0.6 and 13 GeV < i' , 

with 8 data points each, and a reasonable overlap of Zh- 
and i^-distributions. The partial overlap in the z^-cuts is 
used to improve the statistics of the 2 data sets. The full 
data set defined in Section III is characterized by an av- 
erage r intermediate between the 2 above data sets. The 
fitted A for tt-I- production on Ne are plotted in Fig. 5. 
Preliminary HERMES data on a Ne target [4] is com- 
pared with the absorption model of Ref. [13, 14] and the 
energy loss model of Ref. [13, 19]. Though the error bars 
in the 2 subsets are relatively large, experimental data 
hint at a decrease of A with r, which confirms the in- 
terpretation of r as prehadron formation time. Its mod- 
est slope indicates that induced partonic energy loss in 



cold nuclear matter is rather weak, as also predicted in 
Rcf. [16]. 



V. SUMMARY AND CONCLUSIONS 

In this work, I proposed a scaling analysis of hadron 
attenuation in nDIS as a tool to investigate quark 
hadronization in cold nuclear matter, and to distinguish 
parton energy loss from nuclear absorption effects in 
experimental data. The scaling properties of experi- 
mental data and theory computations of the hadron 
attenuation ratio Rm can be summarized by the value 
of the exponent A in the scaling variable r introduced 
in Eq. (3). The exponent A is able to clearly distinguish 
models based on parton energy loss (A < 0) from models 
based on hadron absorption (A > 0.5). Experimental 
data on pion and charged hadron production have 
been shown to scale with r and exhibit A > 0.4, which 
is a clear indication that the hadronization process 
starts on a time scale of the order of a few Fermi, and 
that prehadronic nuclear absorption dominates hadron 
quenching in nuclear DIS. The scaling variable r can 
then be interpreted as a measure of the formation time 



of the prehadron, the color neutral precursor of the 
observed hadron. Note that the scaling analysis cannot 
measure the absolute magnitude of the prehadron 
formation length, only its dependence on Zh and ly. A 
more direct detection of in-medium hadronization, and 
a measurement of the overall scale of the prehadron 
formation time, is possible by looking at the hadron 
PT-broadening, as proposed in Ref. [16]. The scaling 
analysis described in this paper will be a useful cross- 
check of this measurement. Establishing a scaling of the 
prehadron formation time with inverse Q^, as predicted, 
e.g., in Ref. [16], will further constrain the hadronization 
mechanism. A dedicated experimental analysis is needed 
to improve the reach and precision of the scaling analysis 
presented in this paper. 
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